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1. Introduction 

Different cosmological observations provide evidence for the existence of Cold Dark Matter 
(CDM). Through measurements of the WMAP satellite in combination with baryonic acoustic 
oscillation and supernova data [[[]] the relic density of dark matter can be constrained very precisely 
to the interval 

Cl CDM h 2 = 0.1126 ± 0.0036 (1.1) 

at la confidence level, where h in units of 100 km s _1 Mpc -1 denotes the present Hubble expansion 
rate. One possible candidate for dark matter is the so called weakly interacting massive particle 
(WIMP). However, the Standard Model of particle physics does not provide such a particle fulfilling 
all relevant constraints. In contrast, the Minimal Super symmetric Standard Model (MSSM) with 
conserved 7?-parity, contains an appropriate dark matter candidate, the lightest neutralino a 
stable WIMP. 

To calculate the dark matter relic density, the time evolution of the dark matter number density n x , 
described by the Boltzmann equation, has to be solved. As we want to consider the case when 
heavier, unstable supersymmetric particles survive in the Universe for sufficient time to affect the 
dark matter relic density, we have to account for all interactions between two surviving particles i 
and j. This can be described by an effective form [[| f|] 



^ = -3Hn x -(a eS v)[n 2 x -(n e x q ) 2 



(1.2) 



As all particles finally decay into the dark matter particle, the total number density can be written as 
n x = Y,i n i with rii being the number density for each particle species. The effective (co-)annihilation 
cross section is given by 



ij n X n X 

where the sum runs over all MSSM particles i and j. The relative velocity of the two interacting 
particles is described by vy. The ratio between their respective number density in thermal equi- 
librium nj q and the number density of the dark matter particle n x at temperature T, is Boltzmann 
suppressed 

ni; — m, 



eq 

4q ~ ex P 
n X 



(1.4) 



Only if another particle is almost mass degenerate with the dark matter particle, the neutralino in 
our case, the effective cross section gets a significant contribution of the corresponding processes. 
Thus annihilation of two neutralinos, co-annihilation of a neutralino with other gauginos, or co- 
annihilation of neutralinos with light sleptons or squarks can contribute in a sizeable manner. 
Having solved the Boltzmann equation numerically, the neutralino relic density can be calculated 
by 

1 til y/ly 

Pcrit 

where m x is the mass and n x the current number density of the neutralino and p cr jt the critical 
density of the Universe. 

Comparing the theoretically predicted dark matter relic density with the experimentally obtained 
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limits allows one to constrain the MSSM parameter space and thus to combine information of cos- 
mological observations with collider searches and precision measurements. 

However, the theoretical prediction of the neutralino relic density suffers from different sources of 
uncertainties, which arise from cosmology as well as from particle physics, for example. From 
the cosmological point of view the uncertainties are connected to the choice of the cosmological 
model [Q] or the definition of the Hubble expansion rate before Big Bang Nucleosynthesis [g]. 
On the particle physics side, uncertainties arise due to the calculation of essential parameters like 
physical masses or couplings of supersymmetric particles. A possible source for these discrepan- 
cies is for example the different treatment of radiative corrections or a different implementation of 
renormalization group equations in spectrum calculators [Q]. 

Another uncertainty, the one we address in the following, concerns the actual precision of the cal- 
culation of the (co-)annihilation cross sections. Current public dark matter tools like Dark SUSY 
[Q] or micrOMEGAs [Js[] evaluate the relic density on the basis of (co-)annihilation cross sections 
calculated at an effective tree-level. Especially with the even more precise data which will be pro- 
vided by PLANCK in the very near future, the accuracy of the theoretical prediction has to be 
improved to meet the experimental precision. The non-negligible impact of next-to-leading order 
corrections on the dark matter relic density has been discussed in several previous analyses. SUSY- 
QCD corrections to neutralino pair annihilation into a quark-antiquark pair have been studied in 
Refs. [^-[ll]]. The corresponding electroweak corrections have been addressed in Refs. [|l2|-[n|], 
the co-annihilation of a neutralino with another gaugino has been also discussed in Refs. [13, p| ]. 
Further studies regarding neutralino pair annihilation or co-annihilation with a tau slepton, have 
been performed on the basis of an effective coupling approach in Refs. [|l^, 16]. SUSY-QCD cor- 
rections to co-annihilation of a neutralino with a stop have been only considered in Ref. [17] so 
far. Although studying the very specific cases of co-annihilation of a bino-like neutralino with a 
right-handed stop into a top quark and a gluon or into a bottom quark and a W-boson, a significant 
impact on the dark matter relic density was shown. Similar results were found in the aforemen- 
tioned studies. 

However, depending on the region of parameter space, other final states, including those with other 
electroweak gauge and Higgs bosons, can become dominant. Thus we extend the analysis of QCD 
and SUSY-QCD corrections to co-annihilation of a neutralino with a stop by computing the general 
case of neutralino-stop co-annihilation into a quark and a Higgs or an electroweak vector boson. 
In the following, the phenomenology of neutralino-stop co-annihilation in the MSSM is discussed. 
Thereafter the calculation of the next-to-leading order corrections to the relevant processes is de- 
scribed and their impact on the co-annihilation cross sections and the dark matter relic density is 
shown. 



2. Phenomenology of neutralino-stop co-annihilation 

2.1 Neutralino-stop co-annihilation in the context of a 126 GeV Higgs boson 

In certain regions of the MSSM parameter space co-annihilation of the lightest neutralino 
with the next-to-lightest supersymmetric particle (NLSP) can become dominant in comparison to 
other (co-)annihilation processes. A particularly interesting example of such a NLSP is the stop: 
Especially when the trilinear coupling A t gets large absolute values, its chirality eigenstates can 
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Figure 1: Leading-order Feynman diagrams for neutralino-squark co- annihilation into a quark and a Higgs 
boson (0 =/i ,//°,A°,i/ ± ) or an electroweak gauge boson (V — y,7^ ,W^). The w-channel is absent for a 
photon in the final state. 



mix significantly. The lighter mass eigenstate can then be almost mass-degenerate with the lightest 



neutralino, which induces the studied neutralino-stop co-annihilation [18] 



Interpreting the recent observation of a boson with a mass of about 126 GeV [ |22| , |23| ] as a light CP- 
even Higgs boson (h ), a specific choice of parameters in the stop and sbottom sector within the 



MSSM is implied [25]. The reason lies in the interplay of these parameters regarding the lightest 
Higgs boson mass. The leading contribution to its mass arises from a loop containing stops, which 



can together with the tree-level be expressed as follows [|26J, [270 
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(2.1) 



where X t = A t — pL / tan j8 and Msusy = ^J m t x m t 2 ■ The maximal contribution from stop mixing is 
obtained for \X t \ ~ \/6Msusy, which favors a sizable trilinear coupling A t and consequently a rather 
light stop. 

Besides, a light third generation of sfermions is helpful in order to reduce fine-tuning and to stay 
compatible with experimental constraints at the same time, like for example in "natural" SUSY 



models Q20| , [21fl . In the following, we focus on neutralino-stop co-annihilation with a quark and 
an electroweak vector or Higgs boson in the final state. The corresponding leading order Feynman 
diagrams are shown in Fig. |[ At tree level, the process is mediated either by an ^-channel quark, a 
^-channel squark, or a w-channel neutralino or chargino exchange. 



2.2 Parameter study within the pMSSM 

In order to quantify the relative importance of the processes in Fig. [IJ we have performed 
a random scan in the phenomenological MSSM. For simplification we restricted the model to 
the following eight parameters: The slepton sector is characterized by one single mass parame- 
ter M~ t , the squark sector is described by M^, 2 for the first and second generation squarks, and 
the common mass parameter M^ 3 for the third generation squarks. Except for A t in the stop sec- 
tor, all trilinear couplings are set to zero. In our study we parametrize the trilinear coupling as 
T t = Y t A t with Y t being the top Yukawa coupling. With the wino and gluino masses fixed by 
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2M\ = M2 = M3/3, which is motivated by gaugino mass unification at the GUT scale, all gaug- 
ino masses are defined through the bino mass parameter Mi . The Higgs sector is defined by the 
higgsino mass parameter 11, the pole mass of the pseudoscalar Higgs boson via and the ratio tan /3 
of the two vacuum expectation values of the Higgs doublets. According to the SPA-convention 



[ 28 ] these soft-breaking parameters are defined at the scale 2=1 TeV. We have randomly gen- 



erated 1.2 million parameter sets with the eight input parameters lying in the following ranges: 

500 GeV < Mq L2 < 4000 GeV, 200 GeV < M x < 1000 GeV, 

100 GeV <M^ < 2500 GeV, 100 GeV < m A < 2000 GeV, (2.2) 

500GeV<M ( - < 4000 GeV, \ll\ < 3000 GeV, 

\T t \ < 5000 GeV, 2<tan/3<50. 



Using SPheno [29] (version 3.2.1) the mass spectrum and mixing matrices have been deter- 



mined for each point in the parameter space. The contributions of the individual (co-)annihilation 
channels as well as the neutralino relic density Q.%h 2 have been evaluated using micrOMEGAs ^ 
(version 2.4.1). The necessary numerical values of the Standard Model parameters have been 



taken from Ref. [|30p. 

In Fig. |H the relative contribution of the different neutralino-stop co-annihilation channels for the 
generated parameter points as a function of the most relevant input parameters is shown. As can be 
seen, in a large number of parameters sets co-annihilation of a neutralino with a stop contributes 
significantly to the overall (co-)annihilation cross section. 

In order to focus only on experimentally viable scenarios, we impose the following selection cuts 
on the data: 

0.0946 < D, x h 2 < 0.1306, 
120 GeV < m h o < 130 GeV, (2.3) 
2.77 -10~ 4 < BR(b^sy) <4.33-10~ 4 . 

The first selection cut chooses the parameter sets where the relic density lies within a 5 a confi- 
dence interval. Considering a theoretical uncertainty of about 3 GeV, the second one corresponds 
to a medium cut on the lightest Higgs boson mass. The third selection cut is chosen to be a 3a 
interval around the observed value of BR(b — > sy) = (3.55 ± 0.26) • 10~ 4 [pip. The relative contri- 
butions of the studied co-annihilation channels of the remaining parameter sets after imposing the 
cuts are shown in the lower part of Fig. ^[ 

Comparing the data points before and after applying the cuts, several interesting features can be 
observed. The shape of the distribution is not significantly changed, but the density of the points 
is reduced. The main characteristics of the co-annihilation, the degenerate masses of the lightest 
neutralino and the stop, is reflected in the left and left-center column of Fig. ^[ where the depen- 
dence of the relative contribution on the the gaugino mass parameter M\ and the third-generation 
squark mass parameter Mq 3 is depicted. For large values of both input parameters, co-annihilations 
cease to be important and annihilations of stops take their place as the dominating contribution of 
the total cross section. 

After imposing the selection cuts, particularly the channel with the top and the gluon in the final 
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Figure 2: Relative contributions of the neutralino-stop co-annihilation channels for the generated parameter 
points as a function of the input parameters M\, M< l3 , T t , and tan/3 before (top) and after (bottom) applying 
the selection cuts of Eq. (2.4). Shown are the contributions from th° (red), tg (green), tZP (blue), tH° 
(yellow), bW + (cyan), tA° (brown), bH + (pink), and ty (gray) final states. The parameters Mi, Mq 3 , and T t 
are given in GeV. 



state is reduced, whereas the final state with the lightest Higgs remains the dominant contributing 
channel. The reason for this behavior becomes visible in the right-center column of Fig. || where 
the dependence on the trilinear coupling T t is depicted. To match the experimental observations, 
a sizable trilinear coupling T t is preferred, which agrees with our previous discussion of Eq. (|2~T|). 
Moreover, positive values for T t are slightly favored, since they allow a better maximization of the 
Higgs boson mass [25]. This results also in an enhancement of the Higgs-squark-squark coupling, 
which is present in the f-channel of the th° final state and leads to the dominating relative contri- 
bution of the co-annihilation into Higgs final states. This means that the same mechanism which 
drives the lightest Higgs boson mass through important stop-loop contributions to the observed 
value, is also responsible for the increase of neutralino-stop co-annihilation into the lightest Higgs 
boson together with a top quark. 

A similar behavior is found for the Higgs parameter fx, which is not displayed: Larger values of \l 
enhances the Higgs-sfermion-sfermion coupling (mainly for the heavy CP-even Higgs), which has 
the same consequences as for large values of T t . 
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M\ 


%,2 






T, 


m A 


I 1 


tan/3 


m-a 
X, 


m h 


m h o 


m H o 


306.9 


2037.7 


709.7 


1499.3 


1806.5 


1495.6 


2616.1 


9.0 


307.1 


350.0 


124.43 


1530.72 



Table 1: Example for a characteristic scenario within the pMSSM. Given are the input parameters as de- 
scribed in the text and selected particle masses. All values except tan/8 are given in GeV. 



a x h 2 


%% -> th° 


Zfii -> bW + 


x"h -> 'Z° 


Sum 


0.114 


38.5% 


5.9% 


3.4% 


47.8% 



Table 2: Neutralino relic density and relative contributions of neutralino-stop co-annihilation into a quark 
and a Higgs or electroweak gauge boson for the example point of Tab. Ill 



In contrast, the dependence on tan /3 is less pronounced. The remaining input parameters, such as 
those related to first and second generation squarks, sbottoms, and sleptons, as well as the higgsi- 
nos are less important in this context and are therefore not displayed in Fig. ||. 
In the following we present an example for a characteristic scenario with dominant contribution of 
the th° final state. The corresponding input parameters and masses are given in Tab. [j] The value 
of the neutralino relic density and the relative contributions of the dominating channels are listed in 
Tab. ||. For a more detailed discussion of further scenarios we refer the reader to our recent paper 

& 

The chosen scenario is characterized by the dominant neutralino-stop co-annihilation into a top 
quark and a light Higgs boson. This is due to the large value of T t in this parameter point, which 
enhances the Higgs-sfermion-sfermion coupling and thus the squark exchange in the ^-channel. 
Final states including a top quark and a Z-boson as well as a bottom quark and a W -boson con- 
tribute to a clearly lesser extent. Summing up the contributions of all final states, neutralino-stop 
co-annihilation accounts to almost half of the total (co-)annihilation cross section at this point. This 
makes the chosen example point an interesting scenario for the following further studies. 



3. One-loop cross section 

3.1 Corrections and divergence treatment 

The virtual corrections for the co-annihilation processes, shown in Fig. |IJ contain propagator 
corrections (see Fig. vertex corrections (see Fig. |]) and box contributions (see Fig. ||). They 
have been calculated analytically and cross-checked using the publicly available tools FeynArts 
[33], FeynCalc [[34]], and FORM [35]. The arising divergences are regularized by calculating 
in D = 4 — 2s dimensions. To preserve supersymmetry the dimensional reduction regularization 
scheme (DR) has been applied. For canceling the ultraviolet (UV) singularities, corresponding 
counterterms to the relevant MSSM parameters and fields have been introduced. With the choice 
of a hybrid on-shell/DR renormalization scheme, we minimize the sources for potential problems 
connected to sensitive parameters, e.g. the bottom trilinear coupling A/,. Using this scheme we 
can treat consistently the relevant parameters in the quark and squark sector of the MSSM for all 
(co-)annihilation processes over a large region of parameter space. The resulting expressions for 
the virtual corrections and their counterterms have been implemented in a numerical Fortran 
code p6j]. We have explicitly verified that after renormalization all UV divergences cancel. 
A special treatment is necessary for canceling the infrared (IR) divergences, which arise in the 
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qt 




<? 



Figure 3: Self-energy corrections for quarks and squarks at one-loop level in QCD contributing to 
neutralino-squark co-annihilation. 




V/4> £ V/t> V 



flN 



Figure 4: Vertex corrections at one-loop level contributing to neutralino-squark co-annihilation into quark 
and Higgs (0) or electroweak vector (V) boson. The diagram involving the V — g — q ~ q vertex is only 
present for the case of a vector boson in the final state. 




Figure 5: Four-point diagrams at one-loop level contributing to neutralino-squark co-annihilation into quark 
and Higgs (0) or electroweak vector (V) boson. The last diagram is absent for a scalar in the final state. 



diagrams where a gluon is exchanged (see Figs. |] |5]). These divergences cancel against similar 
ones that come from the real radiation diagrams (see Fig. ^ where a gluon is emitted off a quark 
or squark. In contrast to the UV divergences, the IR divergence treatment is slightly more difficult, 
as the divergences of the real emission diagrams come from the integration over the gluon phase- 
space. 

Different approaches exist to deal properly with this kind of divergences. For our study we use the 



phase-space slicing method [{37|-|39p which uses a lower cut on the gluon energy AE in the phase- 
space integration to get finite real corrections. 1 The divergent part of the phase-space integral 
can then be performed analytically in the limit of gluon with low energy - the so called soft- 
gluon approximation. Divergences obtained in the soft-gluon approximation cancel analytically 
exactly with those coming from the virtual collections. The dependence on this introduced cut 



'The implementation of a dedicated dipole subtraction method a la Catani-Seymour | |4c| ] is work in progress and 
subject to a later publication. 
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Figure 6: Real gluon emission diagrams at one-loop level contributing to neutralino-squark co-annihilation 
into quark and Higgs (0) or electroweak vector (V) boson. The last diagram involving the four-vertex is 
absent for a scalar in the final state. 



should in principal completely vanish. However, in practice the cancellation is limited by the 
stability of numerical integration of the real corrections. We have verified that in our calculation 
and Fortran-implementation all cross sections are insensitive to the exact choice of this cut. For 
a more detailed discussion, we refer the reader to our paper [32]. 



3.2 Numerical results 

First, we focus on the impact of the calculated next-to-leading order corrections on the neutralino- 
stop co-annihilation cross section. In Fig. the next-to-leading collection to the cross section av 
(without the tree-level contribution) with its contributions of the different corrections is depicted. 
Although all renormalized contributions are UV finite, the box, vertex and real part of the correc- 
tion are IR divergent. The IR divergent parts still contain uncancelled poles along with uncancelled 
logarithms of the large separation scale, which leads to an ambiguity in the exact definition of these 
contributions. The large remaining logarithms cause the box contribution to be artificially large 
and leads to negative real corrections at the same time. The latter consists of the soft and the hard 
part and is thus cut-off independent. 

In addition one notices that the propagator corrections is significantly enhanced, which holds rela- 
tively also for the box contributions. This is due to the fact that this studied process is dominated 
by the ^-channel exchange. One of the corrections to the ^-channel leads to a correction to the stop 
propagator as well as to a box diagram where a gluon is exchanged between the initial state squark 
and the final state quark. Thus the enhanced box and propagator corrections entail a large overall 
next-to-leading order correction in the case of neutralino-stop co-annihilation into a Higgs boson. 
This is depicted in Fig. [7], where a comparison of our tree-level calculation, the effective tree-level 
calculation of micrOMEGAs and our full one-loop calculation is shown. In the upper panel the 
corresponding cross section av is depicted, in the lower part we show the ratio between the differ- 
ent cross sections. As it can be seen, our tree-level calculation is in agreement with the values given 
by micrOMEGAs. Taking into account our full one-loop calculation the cross section increases by 
about 30% in comparision to the tree-level. This is caused by the significant contribution from the 
box diagrams and propagator corrections as discussed above. 
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Figure 7: On the left-hand side the contribution of the different corrections to the total next-to-leading 
order correction for the process X\h ~^ th° is shown. The gray area indicates the thermal distribution (in 
arbitrary units). On the right-hand side the tree-level (black dashed line), full one-loop (blue solid line) and 
mi cr OMEGAs (orange solid line) cross section for the process X\h — > th° is depicted. The upper part shows 
the absolute value of ov together with the thermal distribution (in arbitrary units), whereas the lower part 
shows the corresponding relative shift (second item in the legend). 

4. Impact on the relic density 

Our numerical implementation of the full on-loop calculation, described in the previous sec- 
tion, can be used as an extension to public dark matter tools like micrOMEGAs to evaluate the 
impact of the one-loop corrections on the neutralino relic density. Even if this study concentrates 
on the case of co-annihilation of the lightest neutralino with the lightest stop, our chosen implemen- 
tation is general so that the code can be used for any neutralino-sfermion co-annihilation process. 
In the following we compare the neutralino relic density obtained by the three in Sec. || discussed 
calculations: the cross-section calculated by default through micrOMEGAs at tree-level, our tree- 
level calculation and our full next-to-leading order calculation. 

First, we study the change of the relic density when varying a single input parameter around the 
value of the studied scenario of Tab. [T[ In Fig. || the relic density Cl x fr is shown as a function of 
the bino mass parameter M\ and the trilinear coupling parameter T t , calculated on the basis of the 
aforementioned three different calculations. 

The predicted relic density is very sensitive to variations of the bino mass parameter. It decreases 
rapidly for higher values of Mi due to a smaller mass splitting between the lightest neutralino 
and the lightest stop, which enhances the stop-stop annihilation along with the neutralino-stop 
co-annihilation. In contrast, slightly lower values for the bino mass parameter enlarge the mass dif- 
ference and suppress the contribution of co-annihilation processes in favor of neutralino-neutralino 
annihilation. The predicted relic density is then higher due to the absence of co-annihilation. 
The lower part of Fig. || shows the relative correction, i.e. the ratio of the relic density calculated 
with our full one-loop co-annihilation cross section to the one included by default in micrOMEGAs 
and our tree-level, respectively, is shown. Taking into account the full one-loop calculation, a 
relative correction to the effective tree-level implemented in micrOMEGAs of about 9% is ob- 
served. This is due to the lightest Higgs final state, which contributes around 38.5% to the total 
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Figure 8: The neutralino relic density £l x h 2 as a function of Mi (left) and T t (right) calculated using different 
co-annihilation cross sections: default micrOMEGAs (orange solid line), tree-level (black dashed line), and 



full one-loop (blue solid line). The gray band indicates the favored range according to Eq. (1.1). The lower 
part shows the relative impact of the one-loop correction on the relic-density compared to the tree-level 
calculation. 



(co-)annihilation cross section with a corresponding correction of around 30% (see Fig. ^). With 
the current experimental uncertainty of about 3% according to Eq. ( |1.1| ), the impact of the presented 
collections is significant and thus important to be taken into account. 

The relic density is less sensitive to varying the trilinear coupling parameter T t around the value of 
the chosen scenario (T t = 1806.5 GeV). This is shown on the right-hand side of Fig. || The dif- 
ference between the uncorrected and corrected relic density in the cosmologically favored region 
corresponds to a difference of approximately 3 GeV in T t . 

More details on the impact of the full next-to-leading order corrections are shown in Fig. ^j, where 
the WMAP favored region is depicted as a function of two parameters. On the left hand side we 
show the dependence on the mass parameter of the third generation of squarks Mq 3 and the bino 
mass parameter M\ . On the right hand side the dependence of the trilinear coupling parameter T, 
and again the bino mass parameter M\ is shown. The WMAP-favored region of parameter space 
within an la interval is marked as a band. In orange we show the cosmologically favored region 
based on the default micrOMEGAs calculation, in blue based on our full one-loop calculation. A 
clear separation of the two bands is visible, which indicates that the resulting corrections are larger 
than the current experimental uncertainties. 

On the left hand side, the relative impact of the one-loop corrections on the relic density is shown 
in black lines. As already discussed before, corrections up to 9% are obtained. Also visible in this 
plot, the WMAP-favored band follows a straight line in the Mi-Mg 3 plane, which corresponds to 
a constant mass difference between the lightest neutralino and the lightest stop of about 40 GeV. 
Above the cosmologically allowed band the neutralino becomes heavier and the mass difference 
decreases. As a consequence, the stop-stop annihilation becomes dominant. As stop-stop annihila- 
tion has typically a significant higher cross section than the co-annihilation processes, this leads to 
a too small neutralino relic density. For large values of M\ (gray area in the upper left corner) the 
stop becomes the lightest supersymmetric particle, which is disfavored as dark matter candidate 
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Figure 9: WMAP-compatible relic density bands from default micrOMEGAs calculation (orange) and our 
one-loop calculation for co-annihilation (blue) in the (Mq 3 ,Mi ) (left) and (If, Mi ) (right) plane. In the plot 
on the left hand side the relative contribution of co-annihilation processes is shown in green contour, and the 
relative impact of the one-loop corrections on the relic density in black lines. The plot on the right hand side 
shows the LSP-NLSP mass difference in green contour, and the lightest neutral Higgs boson mass in black 
lines. 

both for its electric and color charge. Below the cosmologially allowed band, the mass difference 
gets larger and thus the neutralino-stop and stop-stop (co-)annihilation becomes more and more 
Boltzmann suppressed. In contrast, the neutralino annihilation, which has a lower cross section, 
becomes now dominant, which results in too large relic density. 

Focusing on the right hand side of Fig. ||, we see in different green colors the mass difference 
between the lightest and next-to-lightest supersymmetric particle depicted. It confirms that the 
WMAP-favored region follows a contour of a constant mass difference around 40 — 45 GeV. In 
this plot the solid black lines show the mass of the lightest Higgs boson. The whole cosmologically 
favored region lies within the recent Higgs mass limit of 125.2 GeV ±0.9 GeV [|24|]. At the mo- 
ment the cosmologically constraints from WMAP are more stringent than the current bounds on a 
Higgs-like particle. 

We can conclude that the impact of the studied one-loop corrections on the dark matter relic density 
is larger than the current experimental uncertainties by WMAP. Therefore it is necessary to take 
these corrections into account for a solid theoretical prediction of the neutralino relic density. 

5. Conclusions 

A powerful method to constrain the parameter space of theories beyond the Standard Model 
is to compare the predicted dark matter relic density with data from cosmological precision mea- 
surements, in particular from the WMAP satellite. On the particle physics side, one of the main 
uncertainties on the relic density prediction arises from the calculation of the (co-)annihilation cross 
sections of the dark matter particle. 
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Therefore we studied the impact of one-loop corrections to neutralino-stop co-annihilation into 
electroweak vector bosons or Higgs bosons on the neutralino relic density. Especially with the 
recent developments regarding the discoveiy of a new boson with a mass of around 126 GeV, 
which can be interpreted as the lightest Higgs boson within the MSSM, the phenomenology of 
neutralino-stop co-annihilation is even more motivated. A favored sizeable trilinear coupling pa- 
rameter T t increases the relative importance of the neutralino-stop co-annihilation into a top quark 
and a Higgs boson at the same time. In our presented scenario such a process contributes around 
40 % to the overall (co-)annihilation cross section and receives one-loop corrections of around 30 
%, such that an impact on the neutralino relic density of around 9 % is visible. As this is larger than 
the current experimental uncertainty of WMAP, it is necessary to take into account such corrections 
for the theoretical prediction of the dark matter relic density. 

For a more detailed discussion of the impact of neutralino-stop co-annihilation in different scenar- 
ios we refer to our recent publication [32|. The neutralino-stop co-annihilation into a gluon will be 
the subject of a later publication. 
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